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Blazars are a small fraction of all extragalactic sources but, unlike other objects, they are strong 
emitters across the entire electromagnetic spectrum. In this study we have conducted a detailed 
investigation of the broad-band spectral properties of the gamma-ray selected blazars of the Fermi 
LAT Bright AGN Sample (LBAS) . By combining the accurately estimated Fermi gamma-ray spec- 
tra with Swift , radio, NIR-Optical and hard-X/gamma-ray data, collected within three months 
of the LBAS data taking period, we were able to assemble high-quality and quasi-simultaneous 
Spectral Energy Distributions (SED) for 48 LBAS blazars. 



INTRODUCTION 



The Large Area Telescope (LAT) on board of 
the Fermi Gamma Ray Space Telescope, launched 
on 11 June 2008, provides unprecedented sensitiv- 
ity in the 7-ray ban d (20 MeV to over 300 GeV, 



Atwood et al. 2009j | with a large increase oy er its 
predeces sors EGRET iThompson et al. 1993| , and 
AGILE [Tavani et al. 20081 ]. The first three months 
of operations in sky-survey mode led to the compi- 
lation of a list of 205 7-ra y sources with stati stical 
significance larger than lOcr [Abdoetal. 2009aj . 



As largely expected from the previous results, most 
of the hi gh Galactic latitude sources in this catalog are 
blazars [Abdo et al. 2009b] ]. One of the most effec- 
tive ways of studying the physical properties of blazars 
is through the use of multi-frequency data. This ap- 
proach has been followed assembling the Spectral En- 
ergy Distributions (SEDs) of many radio, X-ray, and 
gamma-ray selected blazars. In all cases, however, the 
effectiveness of the method was limited by the avail- 
ability of only sparse, often non-simultaneous, flux 
measurements covering a limited portion of the elec- 
tromagnetic spectrum. 

With Fermi, Swift, and other high-energy astro- 
physics satellites simultaneously on orbit, comple- 
mented by other space and ground-based observato- 
ries, it is now possible to assemble high-quality data 
to build simultaneous and well sampled SEDs of large 
and unbiased samples of AGN. 

We derive the detailed SED of a subsample 
of 48 Fermi blazars using simultaneous or quasi- 
simultaneous data obtained from Swift and other 
ground and space based observatories. 



II. MULTI FREQUENCY OBSERVATIONS 

We describe the multi-frequency observations of 
LBAS blazars carried out between August and Oc- 
tober 2008 with Fermi, and between May 2008 and 
January 2009 with Swift and other space and ground- 
based facilities. 



A. Fermi LAT data analysis and gamma-ray 
energy spectra 

The Fermi-LAT data from 4 August to 31 Octo- 
ber 2008 have been analyzed, selecting for each source 
only pho tons belonging to th e diffuse class (Pass6 
V3 IRF) [Atwood et al. 20091 ]. Events within a 15° 
Region of Interest (Rol) centered around the source 
have been selected. In order to discard photons from 
the Earth albedo, events with zenith angles larger 
than 105° with resp ect to the Earth reference frame 
Abd o et al. 2009a] have been excluded from the data 
samples. 

A maximum likelihood analysis (<7i/ifce)|l|] has been 
used to reconstruct the source energy spectrum. A 
model is assumed for the source spectrum as well 
as for the diffuse background components, depend- 
ing on a set of free parameters. The Galactic diffuse 
emission is modeled using GALPROP package while 
the extragalactic one is described by a simple power 



law Abdo et al. 2009a|. The method has been im- 



plemented to estimate the parameters in each indi- 
vidual energy bin (2 bins per decade, starting from 
100 MeV), and the parameters obtained from the fits 
are used to evaluate the sources fluxes. For each en- 
ergy bin the source under investigation and all nearby 
sources in the Rol are described by one parameter 
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representing the integral flux in that energy bin. The 
diffuse background components are modeled with one 
single parameter each, describing the normalization. 
For each bin, only fit results with a significance larger 
than 3(7 have been retained. Depending on the flux 
and energy spectrum, 4 to 7 bins had positive detec- 
tions for each AGN in the sample. 

Once the differential flux in each energy bin 4>(E) 
has been evaluated, the corresponding SED is then 
obtained by multiplying the differential flux by the 
square of the central energy value of that bin, i.e., 
vF(y) = E 2 cj)(E) where E — hv. The vertical error 
bars represent only the statistical errors. The system- 
atic uncertainties in the effective area for the Pass6 
V3 DIFFUSE event selection have been estimated to 
be 10% at 100 MeV, 5% at 562 MeV and 20% for 
energies greater than 10 GeV [Abdo et al. 2009c| . 



1. UVOT data analysis 

Swift observations are normally carried out so 
that UVOT produces a series of images in each of 
the lenticular filters (V, B, U, UVW1, UVM2, and 
UVW2). 

Counts were extracted from an aperture of 5" , ra- 
dius for all filters an d converted to flux es using the 
standard zero points Poole et al. 200^ |. 

The fluxes were then de-reddened using the appro- 
priate values of E(B — V) for each source taken from 
[Schleeel et al. 1998j | with A X /E(B - V) ratios cal- 
culated for UVOT filt ers using the mea n interstellar 
extinction curve from Fitzpatrick 19 99] . 



2. XRT data analysis 



Swift data analysis 



The Swift Gamm a-Ray-Burst (GRB) Explorer 



Gehrels et al. 2004j is a multi-frequency, rapid re- 



sponse space observatory that was launched on 20 
November 2004. To fulfill its purposes Swift car- 
ries three instru ments on board: the Burst Alert 
Telescope [BAT, iBarthelmv et al. 20"05T ] sensitive in 
the 15-150 keV ban d, the X-Ray Telescope [XRT, 
iBurrows et al. 2005j | sensitive in the 0.3-10.0 keV 
band, a nd the UV and Optica l Telescope (170-600 nm) 
[UVOT. lRoming et al. 20051 ]. The very wide spectral 
range covered by these three instruments is of crucial 
importance for blazar issues as it covers where the 
transition between the synchrotron and inverse Comp- 
ton emission usually occurs. 

The primary objective of the Swift scientific pro- 
gram is the discovery and rapid follow up of GRBs. 
However, there are periods when Swift is not en- 
gaged with GRB observations and the observatory 
can be used for different scientific purposes. The 
sources observed through this secondary science pro- 
gram are usually called Swift fill-in targets. Since the 
beginning of its activities Swift has observed hun- 
dreds of blazars as part of the fill-in program [e.g., 
iGiommi et al. 20071 ]. With the launch of AGILE 
and Fermi, the rate of Swift blazar observations in- 
creased significantly, leading to the observation (and 
detection) of all but 6 blazars in the LBAS sample. 

The Swift database currently includes 119 obser- 
vations of 48 LBAS blazars that were carried out ei- 
ther simultaneously or within three months of the 
Fermi LBAS data taking period. We used the UVOT, 
XRT and BAT data of these observations to build our 
SEDs. Some blazars were observed several times in 
the period that we consider in this paper; in such cases 
we considered only the exposures where the source was 
detected at minimum and maximum intensity by the 
XRT instrument. 



The XRT is usually operated in Auto State mode 
which automatically adjusts the readout mode of the 
CCD detector to the source brightness, in an attempt 
to avoi d pile-up (for details of the XRT observing 
modes) [Burrows" et al. 2005| . Given the low count 
rate of our blazars most of the data were collected us- 
ing the most sensitive Photon Counting (PC) mode 
while Windowed Timing (WT) mode was used for 
bright sources with shorter exposures. 

The XRT data were processed with the XRTDAS 
software package (v. 2.4.1) developed at the ASI Sci- 
ence Data Center (ASDC) and distributed by the 
NASA High Energy Astrophysics Archive Research 
Center (HEASARC) within the HEASoft package 
(v. 6.6.1). Event files were calibrated and cleaned with 
standard filtering criteria with the xrtpipeline task us- 
ing the latest calibration files available in the Swift 
CALDB. Events in the energy range 0.3-10 keV with 
grades 0-12 (PC mode) and 0-2 (WT mode) were used 
for the analysis. 

Events for the spectral analysis were selected 
within a circle of 20 pixel (~47 arcsec) radius, 
which encloses abou t 90% of the PSF at 1.5 keV 
[Moretti et al. 20051 ] , centered on the source position. 
For PC mode data, when the source count rate is 
above ~ 0.5 counts s _1 data are significantly affected 
by pile-up in the inner part of the Point Spread Func- 
tion (PSF). For such cases, after comparing the ob- 
ser ved PSF profile with the analytical model derived 
by [Moretti et al. 2005 1. we removed pile-up effects 



by excluding events detected within up to 6 pixels 
from the source position, and used an outer radius of 
30 pixels. The value of the inner radius was evaluated 
individually for each observation affected by pile-up, 
depending on the observed count-rate. 

Source spectra were binned to ensure a minimum of 
20 counts per bin to utilize the % 2 minimization fitting 
technique. 

We fitted the spectra adopting an absorbed power 
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law model with photon index T x . When devia- 
tions from a single power law model were found, we 
adopted a log-p arabolic law of the form F(E) = 
KE (-a+b-Log{E)) [Massaro et al. 2004j which has been 
show n to fit well the X - ray spectrum of blazars 
[e.g., iGiommi et al. 20051 iTramacere 2009} . This 
spectral model is described by only two parame- 
ters: a, the photon index at 1 keV, and b, the 
curvature of the parabola. For both models the 
amount of hydrogen-equivalent column density (N# ) 
was fixed to the Galac tic value along the line of sight 
[Kalberla et al. 2005| . 



3. BAT hard X-ray data analysis 

We used survey data from the Burst Alert Tele- 
scope (BAT) on board Swift to produce 15-200 keV 
spectra of the blazars presented in this analysis. In 
ord er to do so, we us ed three years of survey data 
see lAiello et al. 20091 for details] and extracted the 
spectra of those blazars that are significantly detected 
in the 15-55 keV band. Because of the very long in- 
tegration time these data are not simultaneous with 
our Fermi data. 

Only 15 blazars, among those presented here, 
were detected by BAT at a significance > 4<r. 
T he spectral extract ion is performed as described 



Aicll o et al. 2008} and the background-subtracted 



spectra represent the average emissions of the sources 
within the time spanned by the BAT survey. 



Swift observations of LBAS blazars carried out 
before May 2008 or after January 2009 



data with other multi-frequency flux measurements 
obtained from a number of on-going programs from 
ground and space-based observatories. 



radio observations 



Quasi-simultaneous radio data for 25 sources of 
the first Fermi bright source catalog were obtained 
within a i^ermz-related monthly broad-band moni- 
toring program including the Effelsberg 100-m ra- 
dio telescope of the MPIfR [F-GAMMA project 
iFuhrmann et al. 20071 1 Angelakis et al. 2008j |. From 
this program, radio spectra covering the frequency 
range 2.6 to 42 GHz were selected to be within the 
time period 4 August 2008 to 31 October 31 2008, 
i.e., quasi simultaneous to the Fermi and Swift ob- 
servations. 



2. OVRO radio data 

Quasi-simultaneous 15 GHz observations of 24 
Fermi LBAS sources were made using the Owens 
Valley Radio Observatory (OVRO) 40 m telescope. 
These observations were made as part of an ongo- 
ing Fermi-LAT blazar monitoring program. In this 
program, all 1158 CGRaBS blazars north of declina- 
tion —20degrees have been observed approximately 
twice per week or m ore frequently since June 2007 
Healev et al. 2008j . For each source, the maximum 



and minimum observed 15 GHz flux densities during 
the 4 August to 31 October, 2008 period were included 
in the quasi-simultaneous SEDs. 



The Swift database includes a number of observa- 
tions of LBAS blazars that were carried out outside 
the period that we consider useful to build our quasi- 
simultaneous SEDs. These measurements are par- 
ticularly important for the case of blazars that have 
never previously been observed by any X-ray astron- 
omy satellite and were below the detection threshold 
of the ROSAT all sky survey. When these Swift ob- 
servations have been analyzed and published by other 
authors we use the flux intensities reported in the lit- 
erature, with particular reference to the latest on-line 
version of the BZcat catalog Q ■ For the cases where 
the Swift results have not yet appeared in the litera- 
ture we estimates the X-ray fluxes from the standard 
pipeline processing that is run at ASDC on all Swift 
XRT data shortly after they are added to the archive. 



3. RATAN-600 1-22 GHz radio observations 

Among the 48 objects for which we present 
Swift and Fermi simultaneous SEDs 32 were ob- 
served between September 10 and October 3, 2008 
with the 600-meter ring ra dio telescope RATAN-600 
[Korolkov &: Pariiskii 1979j of the Special Astrophys- 
ical Observatory, Russian Academy of Sciences, lo- 
cated in Zelenchukskaya, Russia. These observa- 
tions, which produced 1-22 GHz instantaneous ra- 
dio spectra, are par t of a long-term program [e.g., 
iKovalev et al. 2002] to monitor continuum spectra of 
active galactic nuclei with a strong parsec-scale com- 
ponent of radio emission. 



C. Other multi- frequency data 



4- Radio, mm, NIR and optical data from the 
GASP-WEBT collaboration 



In order to improve the quality of our SEDs we com- 
plemented the Fermi and Swift quasi-simultaneous 



In the period considered in this work, the GLAST- 
AGILE Support Program (GASP) originated from the 
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Whole Earth Blazar Telescope (WEBT) [see e.g., 
IVillata et al. 20071 . iRaiteri et al. 2008al | ) carried out 
~ 3000 optical (R band) observations of 19 LBAS 
blazars, while ~ 700 near-IR (JHK, Campo Impera- 
tore), and ~ 600 microwave (230 and 345 GHz, SMA) 
and radio data (5 to 43 GHz, Medicina, Noto, UM- 
RAO) observations were taken on the same sources. 
In the SED plots we report the average, maximum 
and minimum values at each observed frequency in 
the period 4 August - 31 October, 20080]. 



Mid-infrared VISIR observations 



The MIR observations were carried out from 



2006 to 2008 using VISIR jLagage et al. 20041 ] . the 
ESO/VLT mid-infrared imager and spectrograph, 
composed of an imager and a long-slit spectrome- 
ter covering several filters in N and Q bands and 
mounted on Unit 3 of the VLT (Melipal). We per- 
formed broad-band photometry in 3 filters, PAH1 
(A=8.59±0.42 H, PAH2 (A=11.25±0.59 ^m), and 
Q2 (A=18.72±0.88 jum) using the small field in all 
bands (19.2" xl9.2" and 0.075" plate scale). 



Non-simultaneous Spitzer Space Telescope 
observations 



The Spitzer Space Telescope is a 0.85-meter class 
telescope launched on 25 August 2003. Spitzer Space 
Telescope obtains images and spectra in the spectral 
range between 3 and 180 micron through three instru- 
ments on board: the InfraRed Array Camera (IRAC), 
which provides images at 3.6, 4.5, 5.8 and 8.0 mi- 
crons, the Multiband Imaging Photometer for Spitzer 
(MIPS), which performs imaging photometry at 24, 70 
and 160 micron, and the InfraRed Spectrograph (IRS) 
which provides spectra over 5-38 microns in low (R ~ 
60-127) and high (R ~ 600) spectral resolution mode. 
The Spitzer Science Archive include MIPS observa- 
tions of 8 sources belonging to the LBAS sample, all 
of them performed earlier than three months from the 
start of LBAS data taking period. 
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FIG. 1: The SED of 0FGL J0538.8-4403 = PKS0537- 
441. (left) and of 0FGL J0712.9+5034 = GB6 J0712+5033 
(right) 



from the mean flux above 100 MeV, obtained with 
a simple power law source model with fixed spectral 
index —2.1. 



III. BLAZAR SED OBSERVATIONAL 
PARAMETERS 

We now estimate some key observational parame- 
ters that characterize the SED of our blazars, namely, 
the peak frequency and peak flux of the synchrotron 
component (i/jf and ^ffcp)), and the peak frequency 
and flux of the inverse Compton part of the SED {y 1 ^ 
and v^F^)). 

As fitting function we used a simple third degree 
polynomial: 



7. AGILE 7 -ray data 



vF tJ 



a ■ v + b ■ v + c • v + d. 



(1) 



The 7-ray data collected by the Gamma-Ray 
Imaging Detector (GRID) on board of AGILE 
Tavani et al. 2008} for energies greater than 100 MeV 



used in this paper (blue star symbols in the SED fig- 
ures) are extracted from the First AGILE Catalog 
of high-confidence sources detected by AGILE during 
the first 12 mo nths of operations, f rom 9 July 2007 
to 30 June 2008 [Pittori et al. 20091 ]. The differential 
AGILE flux values appearing in the SED figures at 
fixed energy point (E=300 MeV) have been rescaled 



In the case of high redshift sources (e.g., 
J1522. 2+3143), we excluded from the fitting proce- 
dure all points in the optical/UV bands that are likely 
to be significantly affected by Lyman-a forest absorp- 
tion. 

We have estimated the peak of the inverse Compton 
power in the SED {vJp C ) and the corresponding peak 
flux (j/ c F(Vp C )) by fitting the X-ray to 7-ray part 
of the SED, which is dominated by inverse Compton 
emission using the polynomial function of Equation [1] 
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FIG. 2: The SED of OFGL J0722. 0+7120 = S50716+714 
(left) and of OFGL J0730.4-1142 = PKS0727-11 (right) 
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FIG. 3: The SED of OFGL J1512.7-0905 = PKS 1510-089 
(left) and of OFGL J1522. 2+3143 = B2 1520+31 (right) 



There are some objects (e.g. J0722. 0+7120 and 
J1221. 7+2814 ) in which the soft X-ray band is still 
dominated by synchrotron radiation, and only the 
Fermi data can be used to constrain the inverse 
Compton component, so the above method is subject 
to large uncertainties. For this reason, in these cases, 
we have used the ASDC SED [| interface to fit the 
simultaneous data points to a SSC model wit h a log- 
parabolic electron spectrum [Tramacere 2009j |. 

The best fit to both the synchrotron and inverse 
Compton components appear as dashed lines in Figs. 
[TJ [21 [3] and |U The details regarding the estima- 
tion of Spectral Energy Distribution observational pa- 
rameters are accurately reported in the S. Raino et 
al. "Interpretation of blazar SEDs based on broda- 
band quasi-simultaneous observations" Proceedings. 



simultaneous SEDs of a sizable and representative 
fraction of an homogeneous sample (~ 45%) of blazars 
detected during an all sky survey .This collection 
of high-quality, well sampled, nearly simultaneous, 
broad-band SEDs for a large number of blazars is 
unprecedented and allowed us to estimate a number 
of important parameters characterizing the SED of 
selected blazars and to address some key aspects of 
blazar demographics and physics. Indeed we con- 
clude that the distribution of the synchrotron peak 
frequency is very different for the FSRQ and BL Lac 
subsamples. In fact for FSRQs it starts at ~10 12 ' 5 
Hz, peaks at ^10 133 Hz and it does not extend be- 
yond k 10 14 5 Hz, the distribution of BL Lacs is much 
flatter, starts at ~10 13 Hz and reaches much higher 
frequencies ( « 10 17 Hz) than that of FSRQs. 



IV. CONCLUSIONS 

We have carried out a detailed investigation of the 
broad-band (radio to high-energy) spectral properties 
of the LBAS sample of Fermi bright blazars using 
a large number of multi-frequency simultaneous ob- 
servations as well as literature and archival data. Us- 
ing data obtained with Fermi, Swift, radio/mm tele- 
scopes, infra-red, and optical facilities. For the first 
time high-quality quasi-simultaneous SED of blazars 
are available for a considerable number of blazars and 
this subset is representative of the entire LBAS. We 
have been able to assemble simultaneous or quasi- 
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FIG. 4: The SED of 0FGL J1159.2+2912 = 4C29.45 (left) 
and of 0FGL J1221. 7+2814 = ON231= W Comae (right) 
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